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Abstract 


The objective of the present study is to examine the potential of using the near-surface 
total chlorophyll a concentration (C sur f), as it can be derived from ocean color 
observation, to infer the column-integrated and the vertical distribution of the 
phytoplanktonic biomass, both in a quantitative way and in a qualitative way (i.e., in 
terms of co mmunit y structure). Within this context, a large HPLC (High Performance 
Liquid Chromatography) pigment database has been analyzed. It includes 2419 
vertical pigment profiles, all sampled in Case-1 waters with various trophic states. 
The relationships between C su rf and the total chlorophyll a vertical distribution, as 
previously derived by Morel and Berthon [1989], are fully confirmed, as the present 
results coincide with the previous ones. This agreement allows to go further, namely 
to examine the possibility of extracting relationships between C SUI f and the vertical 
composition of the algal assemblages. Th anks to the detailed pigment composition 
available from HPLC measurements, the contribution of three size classes (micro-, 
nano-, and pico-phytoplankton) to the local total chlorophyll a concentration can be 
assessed. Corroborating previous findin gs ( e.g ., large species dominate in eutrophic 
environments, whereas tiny phytoplankton prevail in oligotrophic zones), the results 
lead to a statistically based parameterization. The predictive skill of this 
parameterization is successfully tested on a separate data set. With such a tool, the 
vertical total chlorophyll a profiles associated with each size class can be inferred 
from the sole knowledge of C S urf- By comb inin g this tool with satellite ocean color 
data, it becomes conceivable to quantify on a global scale the phytoplankton biomass 


associated with each of the three size classes. 


1. Introduction 


A permanent and global monitoring of the world ocean, and in particular of the algal 
content of its upper layer, can only be achieved by satellite observation of the “ocean color”. 
Ocean color, as remotely detected, provides the chlorophyll concentration restricted to the 
upper layer only [ Gordon and Mccluney , 1975]. Within the context of ecological studies 
dealing with the vertical distribution of algal species, as well as for bio-geochemical 
applications involving primary production, such satellite information about the upper layer is 
insufficient. Indeed, the assessment of the algal biomass must be extended downward, in 
order to encompass the entire column where algae can live and grow. For instance, the 
transformation of “chlorophyll maps” as obtained from spacebome sensors into “primary 
production maps”, through the use of a light-photosynthesis model, requires at least that the 
column-integrated) vegetal biomass is known, and better, that the biomass vertical distribution 
within the decreasingly illuminated layers can be in some way described. The terms “known” 
or “described” actually mean “assumed” or “predicted” with a sufficient degree of 
confidence, and for each “pixel” of a satellite image. 

The need for extending the upper layer chlorophyll concentration was part of the 
motivation of the first statistical study presented by Morel and Berthon [1989, hereafter 
denoted MB89], which dealt with the shape of the phytoplankton vertical distribution and its 
possible relationship with the near-surface chlorophyll a concentration. This previous study 
was based on the examination of about 4000 vertical profiles of chlorophyll a in Case-1 
waters only, and the main conclusions of this analysis were as follows: 

i) The integrated biomass over the euphotic layer, and the surface concentration, are 
highly correlated, but in a nonlinear fashion. 

ii) The vertical profiles in stratified waters can be sorted according to the surface 
chlorophyll a concentration into several trophic regimes (or trophic categories), and each of 
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in) these situations exhibits a typical shape, which includes a more or less pronounced, 
and more or less deep, maximum; for such stratified waters, a parameterization was proposed 
which allows the chlorophyll a vertical profiles to be predicted in a continuous manner from 
the surface chlorophyll a value. 

iii) And finally, well-mixed waters ( i.e ., when the pycnocline is deeper than the euphotic 
depth) exhibit, as expected, substantially uniform chlorophyll a profiles, and have to be 
considered separately. 

With the same general aims, the present study is firstly motivated by the considerable 
methodological improvement brought by the introduction, and then the systematic use, of the 
High Performance Liquid Chromatography (HPLC) method. Thanks to this unique technique, 
the determination of the chlorophyllous pigment concentration is more accurate than with the 
previous, sometimes diverging, fluorometric and spectrophotometric techniques, as faced 
with m MB89 study. The first and logical goal is to check whether the global description 
given in MB89 remains valid (confirmed or modified) when another data set, based on 
another analytical method and a more comprehensive sampling of the world ocean, is 
considered. 

Beside the point of revisiting the previous parameterizations, which were only 
expressed with respect to the chlorophyll a, it becomes increasingly obvious that the use of 
this single pigment as a proxy of the algal biomass remains insufficient as far as oceanic 
biogeochemical cycles are to be studied and ultimately modeled. In effect, the “quality” of 
the phytoplankton population (namely its taxonomic composition) impacts on, or 
reciprocally, is a signature of, specific biogeochemical processes. For example, tiny 
phytoplankton are preferentially associated with the presence of regenerated forms of the 
nutrients they are able to utilize, whereas large phytoplankters (diatoms), which are more 
involved in so-called new production, develop preferentially when fresh nutrients become 
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available [Eppley and Peterson, 1979; Malone, 1980; Goldman, 1993]. A relationship 
between the trophic status of an oceanic system, as revealed by its near-surface chlorophyll a 
content, and the taxonomic composition of the whole algal assemblage can thus be 
reasonably expected [ Claustre , 1994]. It is, therefore, timely to investigate whether a 
meaningful indication regarding the phytoplanktonic community structure can be inferred 
from the sole knowledge of the chlorophyll a concentration within the near-surface layer. The 
analysis of an HPLC pigment database offers an invaluable possibility in this direction. 

As mentioned earlier, the HPLC method has brought an improved assessment of the 
total chlorophyll a concentration — for the sake of simplicity, the complex vocable “total 
chlorophyll a concentration” will be noted C; see Table 1 for a list of symbols — . In addition, 
and perhaps more importantly, this method allows a suite of accessory pigments (carotenoids 
and chlorophylls) to be determined. Many of these pigments are specific of individual 
phytoplanktonic taxa or groups [. Jeffrey and Vesk, 1997]. They can thus be used as 
biomarkers, and eventually assigned to different phytoplankton size classes, such as micro, 
nano- and pico-pbytoplankton [Vidussi et al., 2001]. Furthermore, the vertically integrated 
biomasses for each algal class, or their vertical profiles, can now be determined via HPLC 
analysis, so that their respective contributions to the total standing algal stock can be 
assessed. The second aim of the present study is, therefore, to examine, from the analysis of 
HPLC data, whether some generic properties regarding the composition and vertical 
distribution of the phytoplankton assemblages may be inferred from the sole total chlorophyll 
a concentration in the near-surface layer, C sur f. 

To achieve these objectives, a large HPLC pigment database has been analyzed. It 
encompasses stations sampled in various trophic and hydrological conditions encountered in 
the open ocean, and all in Case-1 waters. The data include vertical profiles of total 
chlorophyll a as well as of accessory pigments. A statistical analysis of this database is 
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performed to extract, and parameterize the relationships between the trophic status, as 
depicted by C SUI f, and (i) the total chlorophyll a vertical distribution, as done in MB89, and 
(ii) the phytoplanktonic community composition and vertical distribution. 


2. Data 

Two independent pigment databases are used in the present study. Both comprise data 
from a large diversity of cruises to the open ocean exclusively. The first one includes 2419 
vertical profiles, and the second one includes 4238 samples restricted to near-surface waters 
only. The former is used to perform the statistical analysis. The latter is exclusively used for 
validating the statistical relationships between C SUI f and the phytoplankton composition within 
the surface layer. The organization and use of the databases are presented in Figure 1 . 

2.1. Data sources 

The geographic distribution of the stations included in the database #1 is uneven 
(Figure 2). The Mediterranean Sea and the Atlantic and Pacific oceans are well represented 
(especially the subtropical gyres of the northern basins), but the South Pacific Subtropical 
Gyre and the Indian Ocean are inadequately represented. It must be noted that there is no data 
for the Arctic Ocean, and most of the stations in the Southern Ocean actually are located in 
the Ross Sea. 

The geolocations for the second database are principally in the Atlantic Ocean. The 
data were collected either during cruises not belonging to the first database, or, if during the 
same cruises, at locations differing from those having provided the vertical profiles data (e.g., 
surface samples taken along transects between stations). 

A summary of the data collection activities for both databases ’ is presented in 
Appendix A (see electronic supplement). 
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2.2. Overall quality control 

As these samples were collected by several teams, and were analyzed in different 
laboratories by using a variety of instruments and HPLC methods, a quality control is 
required in view of coherently homogenize these data sets. 

The first step of the quality control applies to all individual samples, those from the 
surface and those belonging to vertical profiles. The detection limit of the HPLC method 
depends on the sensitivity of the equipment and on the filtered volumes. Samples where total 
chlorophyll a concentration, C, is below 0.001 mg m' 3 are rejected; for accessory pigments 
concentrations (denoted P for any considered pigment) values below 0.001 mg m’ 3 are reset 
to zero. This rejection has no real impact on the amount of significant data, because such very 
low pigment concentrations are normally encountered at great depths (i.e., beyond the depths 
under investigation). 

As shown by Trees et al. [2000], C and the sum of the concentrations of major 
accessory pigments are tightly correlated, and covary in a quasi-linear manner. The same 
kind of covariation is thus expected, and actually does exist, m the two databases used here. It 
can be used as a tool for a statistical test to identify and eliminate deviant data. In this study, 
the rejection rate is purposefully not severe, as it leads to eliminate the outliers 
(approximately 13%) for which the deviation exceeds three standard deviations with respect 
to the mean. Indeed, a balance must be maintained between the need to remove unwanted 
experimental errors, and the need to maintain enough data to respect the likely natural 
variability of the database. 

The second step of the quality control deals only with the pigment vertical profiles, 
which are checked one by one. The systematic criteria adopted when selecting the profiles are 
as follows: (i) only the profiles where the euphoric depth, Z cu , was reached are kept; (ii) the 
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uppermost sample of the profile must have been collected between the surface and 10 m 
depth; and (iii) each profile must include a minimal number of samples, allowing a proper 
description of the shape of the vertical profiles. Actually a visual inspection of each profile is 
necessary to eliminate those imperfectly described because of a misconceived sampling 
strategy or failures in the experiment. In total, this quality control process has led to the 
rejection of about 1 8% of the initially available profiles. 

2.3. From accessory pigment concentration to phytoplankton size classes 

The HPLC pigment analysis allows the determination of many compounds (up to 15 
different pigments); some of them are redundant in term of their chemotaxinomic 
significance. Thus, it is appropriate to condense the information contained in these 
multivariable data. Along the lines of previous studies [ Claustre , 1994; Vidussi et ai, 2001], 
“pigment indices” or pigment-based identifiers are constructed with the objective of 
extracting the quantitative information (concentration) and the qualitative information 
(taxonomic composition) encapsulated within the full suite of pigments. For this purpose, 
seven pigments are selected because they provide efficient discriminating tools to identify 
particular phytoplanktonic classes among the diversity of the oceanic algal community. These 
seven pigments are: fucoxanthin, peridinin, 19'-hexanoyloxyfucoxanthin, 19'- 

butanoyloxyfiicoxanthin, alloxanthin, chlorophyll b and divinyl-chlorophyll b, and 
zeaxanthin. A list of these seven pigments, along with their taxonomic significance and their 
relationship with algal size are summarized in Table 2. 

A multiple regression is performed between the total chlorophyll a content integrated 
over the euphoric layer, (C)^ , and the content of each of the seven pigments integrated over 
the same layer, (T) Zeu ■ This regression is highly significant (r 2 = 0.76; n = 2419; p < 0.001). 
The coefficients of this regression (Table 2) are estimates of the average ratio of total 
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chlorophyll a to each accessory pigment [see also Gieskes et al., 1988; Everitt et al., 1990; 
Bustillos-Guzman et al . , 1995]. These slopes are used to form the weighted sum of the 
diagnostic pigments concentrations, XDP, expressed as: 

XDP = 1.41 Fuco + 1.41 Peri + 1.27 19'-HF + 0.35 19'-BF + 0.60 Alio 

+ 1.01 TChl6 + 0.86 Zea (1) 

Note that XDP corresponds to the C value which can be reconstructed from the knowledge of 
the concentration of the seven other pigments. Now, the relative proportion (%) of the 
biomasses associated to the three phytoplankton size classes are subsequently derived 
according to 


[%micro] = 100 (1.41 Fuco + 1.41 Peri) / XDP (2a) 

[%nano] = 100 (1.27 19'-HF + 0.35 19'-BF + 0.60 Alio) / XDP (2b) 
[%pico] = 100 (1.01 TChl& + 0.86 Zea) / XDP (2c) 


These three exclusive groupings represent the aforementioned pigments indices. It is 
straightforward to verify that their sum is equal to 100%. The total chlorophyll a 
concentration associated with each of the three size classes is then derived according to: 


micro-C = C [%micro] 

(3a) 

nano-C = C [%nano] 

(3b) 

pico-C = C [%pico] 

(3c) 
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A recent HPLC mtercomparison exercise of pigment determination on natural 
samples, involving four laboratories, showed that C can be determined with an accuracy to 
within 8% [ Claustre et al, 2004]. Furthermore it was shown that accessory pigment ratios, as 
those in equations (2), are determined with a better accuracy than individual accessory 
pigment concentrations (because the normalization processes cancel part of the analytical 
bias). Thus it is expected that the use of equations (3) is very likely the most accurate 
approach when dealing with different data sources. 

Compared to the previous parameterization proposed by Vidussi et al. [2001], the 
present one is improved in that it explicitly accounts for the differing total chlorophyll a to 
accessory pigment ratios, typical of each major phytoplankton group (all set to unity in 
Vidussi et al. [2001]). 

2.4. Using and sorting the data 

All of the 2419 C vertical profiles included in the database #1 (see flowchart in Figure 
2) are first used to study the relationship between C SU rf and (C'j . In a second step, these 

profiles are also used to develop a parameterization of the shape of the C profiles according 
to the trophic status, as depicted by C sur f. This phase of the work tests the MB 89 results 
against an entirely different database. The stratified waters where C sur f is greater than 4 mg m' 
3 are incorporated into the mixed waters category. With such concentrations, mdeed, Ze U is < 
20 m, leading to rather uniform vertical distributions, similar to those found in mixed regime 
(see later). 

In a second phase, the same database #1 is split into two subsets using a random 
process. The first subset comprises 80% of the vertical profiles (1936 stations), to be 
statistically analyzed for deriving the vertical average profiles of fractional C (micro-C, nano- 
C, and pico-C), and proposing a parameterization of these profiles. The remaining 20% (483 
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stations) are kept separately to validate the parameterizations. Actually, a validation is needed 
because such a parameterization of the fractional C profiles has never been attempted, and the 
validation requires a database which includes the full suite of pigments. 


3. Methods 

3.1. Analysis of the total chlorophyll a vertical profiles: methodological background 

The present statistical analysis follows the scheme described in MB89. This method 
was developed for the sum of chlorophyll a plus phaeopigments a (denoted “C” in MB89), 
and will now be used for total chlorophyll a in this study, (also denoted C, as mentioned 
earlier). The same method will be extended without any change to any other pigment or 
group of pigments. Very low concentrations of phaeopigment are usually detected in oceanic 
waters, when determined by HPLC, and can actually be neglected. 

The procedure firstly requires that the stations sampled in well-mixed waters and in 
stratified waters be separately considered. This sorting is made according to the ratio between 
the depth of the euphotic layer, Z eu , and the depth of the mixed layer, Z m , which is taken from 
the World Ocean Atlas 1994 [ Monterey and Levitus, 1997]. When Zeu/Zm is < 1, the sample 
involved is assumed to be from a well-mixed water column, and when Zeu/Zm is > 1, from a 
stratified water column. Z eu was not usually available as an independent measurement, so it 
was inferred from the C vertical profile itself, using a bio-optical model. In MB89, the model 
of Morel [1988] was used. This model has been recently revised [Morel and Maritorena, 
2001], yielding only minor changes (actually slightly larger estimates) to the derived Zeu 
values in oligotrophic waters, and insignificant change for other waters. This new model is 
used in the present study. 
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The statistical analysis is then performed on the basis of two dimensionless quantities. 
A dimensionless depth, £, is introduced, which is obtained by dividing the actual 
(geometrical) depth, z, by Z cu : 


C=z/Z eu (4) 

Even when scaled with respect to the various profiles with differing pigment content 
cannot be straightforwardly compared. Another normalization is thus needed, which consists 
of considering dimensionless concentrations, c(Q. These concentrations are obtained by 
dividing C(Q, the local (in terms of Q concentration, by Czcu, the average concentration 
within the euphotic layer, which is obtained as 


Zcu 

Czcu=(Z cu )- 1 jc(z)dz, (5) 

0 

and 

c(0 = C(Q / Czcu . (6) 

After the aforementioned double normalization has been applied, the dimensionless c(Q 
profiles can be pooled together, and compared regardless of their absolute magnitudes. 

In anticipation of further statistical analyses of the profiles and their shapes, it is 
necessary to make the data easy to handle, i.e., with regularly spaced c(Q values along the 
dimensionless depth. Consequently, for each station, a linear interpolation (with a constant 
increment in set to 0.1) is performed between the data of the actual profiles. This 
interpolation is extended down to twice Z eu (^ = 2), whenever possible. The rationale for this 
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extension lies in the frequent occurrence of deep chlorophyll maxima, often located in 
oligotrophic waters in the vicinity of, or slightly beyond, Z eu [ Letelier et al, 2004], 

The dimensionless c(Q profiles are then sorted according to their absolute C su rf value 
into different “trophic categories”, defined by successive intervals within the C sur f continuum. 
For each category, an average c(Q profile is computed, with its standard deviation (see 
Appendix B in electronic supplement). 

Note that when dealing with any other pigments, or groups of pigments (like pico-C, 
nano-C or micro-C) the same kind of normalization can be performed: the particular local 
pigment (or group of pigments) concentration, P(0, is again normalized with respect to Czcu, 
according to 


P(0 =P(0 / Czcu (6’) 

so that p(Q represents a relative proportion of a given pigment with respect to the average 
concentration of total chlorophyll a within the euphotic zone. 

3.2. Total chlorophyll a vertical profiles parameterization 

The regular evolution of the c(Q profiles as a function of C sur f (discussed below) 
allows its parameterization. In MB89, each average c(Q profile was modeled by using the 
generalized Gaussian profile [Lewis et al., 1983], superimposed to a background 
concentration; thus, a profile was described by four parameters (three for the Gaussian profile 
and one for the background). A slightly modified version of this model is used here, with the 
objective of improving the representation of both the surface and the deepest values, the latter 
being generally lower than the former. This modification consists of replacing the constant 
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background with a linear function that decreases downward from a surface value (Cb). This 
substitution requires the introduction of a fifth parameter to fix the slope (s) of this decrease: 

c(Q = C b - s C + C_ exp {-[(?- ^ ) / AG 2 }, (7) 

where represents the maximum concentration, ^max is the depth at which this maximum 
occurs, and depicts the width of the peak. The use of equation (7) will is explained below 
(Section 4.1.2.). 


4. Results and Discussion 

4.1. Total chlorophyll a biomass 

4.1.1. Column-integrated total chlorophyll a 

In MB89, an initial distinction was made between Cpd, the arithmetic mean 
concentration over the first penetration depth [i.e. Z cu /4.6, Gordon and Mccluney, 1975] and 
C sat , the weighted concentration within the same layer. This distinction has been shown (in 
MB89) to be useless, as these quantities are in effect confounded. They are replaced here by a 
unique quantity (denoted C SUI f, as mentioned earlier) which represents the arithmetic average 
over the first penetration depth. The column total chlorophyll a content integrated over the 
euphotic zone, (C) 2m , is then studied as a function of Csurf- The log-log plots of these 
quantities (Figure 3) suggest nonlinear relationships of the form 


are appropriate (as in MB89), and that the exponent B is notably below unity. 
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In the plot for stratified regimes (Figure 3a), a distinct change in the slope seems to 
occur around C sur f ~ 1 mg m' 3 , so two separate regression analyses on each side of this 
threshold are performed (Table 3). Note that a similar discontinuity was already noted in 
MB89. In the domain C sur f < 1 mg mf 3 , the scatter of the points is more important than in the 
C SU rf > 1 mg m' 3 domain, and reflects the variability associated with the existence of a deep 
chlorophyll maximum typical of oligotrophic conditions. 

In mixed waters, the log-log plot of (C)^ vs Csurf (Figure 3b) shows a unique linear 

trend, a steeper slope, and a low scatter in the data, as a result of the rather homogenous 
pigment distribution. The corresponding power law is highly significant (Table 3). 

The present results are very consistent with, and independently confirm those of 
MB89 (the dashed lines in Figure 3). The exponent in the relationship for stratified waters 
and when C SU rf S 1 mg.m' 3 is slightly lower than in the MB89 relationship (0.357 vs 0.425), 
which leads, in oligotrophic waters, to higher (C)^ values than in the MB89 study. The 
differences are statistically insignificant. For mixed waters, the present relationship is 
essentially identical to that of stratified waters with C sur f > 1 mg mf 3 . The closeness of the 
exponents for stratified waters when Csurf 5 " 1 mg m~ 3 and for mixed waters (0.615 vs 0.538 
respectively) was also reported in MB89. It emphasizes that (notwithstanding the deep 
stratification) the algal biomass in eutrophic conditions remains homogeneously distributed 
within the euphotic layer, essentially as a consequence of its reduced thickness. 

4.1.2. Total chlorophyll a vertical profiles 

On the basis of the C SU rf values, the vertical C profiles have been sorted as in MB 89, 
except that more trophic categories, namely nine for the stratified waters (SI to S9), and five 
for the mixed waters (Ml to M5) have been presently considered, instead of seven and two in 
MB89, respectively. The C sur f intervals selected for defining each category axe provided in 
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Table 4; note also that the lower limits are very different, with < 0.04 mg m' 3 for SI, and < 
0.4 mg m' 3 for Ml. This more detailed categorization scheme was made possible because the 
d a t a used in this study are rather uniformly distributed, which permitted a significant number 
of profiles to be enlisted in each of the nine or five trophic categories. 

Stratified waters 

In the stratified regime, the shapes of the average dimensionless profiles exhibit a 
remarkably regular evolution along the trophic status (Figure 4a). When C sur f increases, the 
change in shape is characterized by i) a smooth ascent of the C maximum, from a depth 
corresponding to about Z cu for SI, up to the near-surface for S9; ii) a regular decrease in its 
magnitude; and iii) a widening and weakening of the maximum that tends to vanish in 
eutrophic waters. It is worth noting that the regular changes of the average profiles are similar 
to those described in MB89 (their categories a to g, for stratified regimes). 

The average dimensionless profiles are used in conjunction with equation (7). The 
fitting procedure (a Newton-type algorithm from the library of the R software [R 
Development Core Team, 2004]) allows the five parameters involved in equation (7) to be 
derived for each category (for the corresponding average C su rf value), from SI to S9 (Table 
4). The modeled profiles (Figure 4b) compare well with the average profiles. The values of 
the five parameters for each of the nine trophic categories are tabulated in Table 5. 

The regular evolution in shape, as described above, suggests the parameterization can 
be used in a continuous fashion. For any Csurf value, the values of the five parameters can be 
interpolated between the discrete values specific of each stratified waters category. This 
interpolation can be computed linearly (present study) or by using fitted polynomials (as done 
in MB 89). 
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The dimensionless profiles can be easily restored to their dimensional concentrations 
and depths, by using the actual Z eu and Czeu values. Both are retrieved using the relationships 
of equation (8), linking (C)^ to C surf , and Z eu to (C)^ [equation (6) in Morel and 
Maritorena, 2001]. 

These restored (or “rescaled”) profiles (Figures 4c, d) provide a direct visual appraisal 
of the evolution in the true profiles when the surface concentration changes. For instance, 
with a Csurf value increasing from about 0.03 to 3.0 mg m~ 3 (SI to S9), the chlorophyll 
maximum increases from 0.17 to 3 mg m' 3 , whereas Z eu decreases from 120 m to 26 m. 

Mixed waters 

In the mixed regime, the average normalized C profiles (Figure 5a) show a rather 
featureless shape: no maximum (except for Ml), a constant concentration which remains 
approximately equal (±10%) to the average concentration within the euphotic layer, and 
finally a systematic decrease of the average concentrations (by about a factor of 2) between 
Z eu and 2 Z cu - The values observed at 2 Z^u, however, are still about 40% of the average value 
within the euphotic column, instead of 20% in stratified waters. The rescaled profiles 
(Figures 5b and c) show that C varies from approximately 0.15 to about 6-7 mg m 3 . The first 
situation (Ml) is encountered during the winter vertical mixing in the North Atlantic, and 
also in the Mediterranean Sea, where the mixed layer depth can reach 200 m [with C « 0.15 
mg m' 3 ; Marty et al., 2002], The category M5 is typical of upwelling conditions, with a 
mixed layer extending deeper than the thin euphotic layer, and C ranging from 4 to about 30 
mg m~ 3 . A remarkable agreement between the present study results and those obtained in 
MB89 is noticed for these mixed regimes. 
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4.1.3. Closer comparison with the MB89 parameterization 

For such a comparison, forty C SUT f values, regularly spaced in terms of their 
logarithms, have been generated within the interval 0.04 and 3 mg m' J . For each of these C sxirf 
values, two C profiles are reconstructed by using the MB89 parameterization and the present 
one. Then, three main characteristics of these profiles are compared: the average 
concentration, (C)^ , the peak value, C max , and the actual depth of this peak, Z max (Figure 6). 
Globally, both parameterizations generate similar values for these shape characteristics. The 
only significant divergence is for Z inax which, in the middle range (20-30m), is smaller with 
the present parameterization than it was by using MB89. This is partly a result of the 
differing bio-optical models used to infer Z eu . For mixed waters (without a deep maximum), 
there is no equivalent comparison; the comparison actually reduces to the relationship 
between (C) Ztu and C sur f, which has already been recorded as being identical in the two 
studies (Figure 3b). 

4.2. Phytoplankton composition and vertical distribution 

Only 80% of the database #1 is now used to study the phytoplankton composition and 
its vertical distribution (Figure 2). A required first step in this analysis is to verify that this 
subset represents the whole set without significant bias. With this objective in mind, the same 
statistical analysis for stratified and mixed waters is repeated on this partial data set and 
compared to the results obtained when using the whole set. Regarding the four characteristics 
representative of each trophic category (i.e., the average C SUi f, Czeu, (C) Zeu , and Z eu ), no 

significant change can be detected between the subset and the entire set — the average values 
and the associated standard deviations for each of the above characteristics, obtained for the 
subset, deviate by less than 1% from those derived for the entire data set. Consequently, the 
following analysis based on the 80% subset can be safely undertaken. 
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4.2.1. Micro-, nano-, and pico-phytoplankton in stratified waters 
Column-integrated contents 

By using the detailed pigment composition together with equations (2) and (3), the 
fractional C, related to each of the three size classes, can be computed at each depth and for 
each profile. These profiles are vertically integrated down to the euphoric depth, and then 
averaged within each trophic category (SI to S9). The comparison of these column- integrated 
fractional C with (C) Zru (all populations confounded) is shown as relative proportions (in 
percent), or as absolute values (in mg m' 2 ), in Figures 7a and 7b respectively, with both 
plotted as a function of C sur f. 

In the most oligotrophic waters (C sur f ~ 0.03 mg m‘ 3 ), pico- and nano-phytoplankton 
are of equal importance; all together they amount to 90% of the total chlorophyll a integrated 
content. In contrast, micro-phytoplankton dominate (75%) in eutrophic waters (Csurf « 3 mg 
m~ 3 ), whereas nano- and pico-phytoplankton do not represent more than 20% and 5% of 
(C) 2 cu ’ respectively. In mesotrophic waters (C su rf « 0.3 mg m" 3 ), nano-phytoplankton are 

predominant (50%), while pico- and micro-plankton account for the remainder in similar 
proportions (25% each). When the absolute contents (mg m" 2 in Figure 7b) are considered, a 
simple picture emerges; the increase in the C values along the trophic continuum (two orders 
of magnitude in C su rf) is essentially caused by the increase in the micro-phytoplankton 
biomass (by a factor of seven); meanwhile the pico-phytoplankton population appears as 
being a rather constant background (» 5 mg m' 2 ), whereas the nano-phytoplankton population 
experiences an increase (actually a tripling) from oligotrophic to eutrophic conditions. 

The observed variations in the phytoplankton assemblages, here quantified as a 
function of the trophic gradient, are in agreement with previous microscopic observations and 
ecological studies. It is widely acknowledged that the phytoplankton size is associated with 


17 


particular types of trophic environment [ Malone , 1980; Chisholm, 1992]; pico-phytoplankton 
would be the main biomass contributor in extremely oligotrophic systems (e.g., subtropical 
gyres), whereas diatoms, which dominate the micro-phytoplankton class, are predominant in 
eutrophic areas (e.g., upwellings). Furthermore the present quantitative analysis reveals the 
rather global ubiquity and stability (20-50% of the biomass) of the nano-phytoplankton 
biomass, which is in agreement with previous observations of 1 9'-hexanoyloxyfucoxanthin 
(the main pigment contributing to nano-C) being an ubiquist pigment in many oceanic 
situations [Ondrusek et al., 1991; Vidussi etal, 2001]. 

When the relative proportions are analyzed for the surface layer only (Figure 7 c) and 
compared to those recorded for the euphotic layer (Figure 7a), noticeable differences are seen 
for the S2-S4 categories: the pico-phytoplankton proportion is enhanced at the expense of the 
nano-phytoplankton. 

Vertical Profiles 

In stratified waters, the average dimensionless profiles of micro-C, nano-C, and pico- 
C exhibit significant and regular evolutions according to the trophic regime (Figures 8a-c). 
The profiles relative to the pico- and nano-phytoplankton classes hold much m common with 
those for the total phytoplankton population (Figure 4a), and are characterized by a distinct 
deep maximum becoming shallower and smoother as C sur f increases. This maximum is 
slightly sharper for pico- than for nano-phytoplankton in most oligotrophic waters (S 1 to S4); 
nevertheless, the peaks have approximately the same intensity. Note that for the most 
oligotrophic category (SI), the dimensionless concentrations of pico- and nano- 
phytoplankton are each about equal to Cz«,, which is possible because the corresponding 
total chlorophyll a value is about twice Cz™ (cf. Figure 4a). Globally, the deep C maximum 
for these waters is almost entirely made up of pico- and ‘nano-phytoplankton assemblages 
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(each with about 40%), in agreement with the column-integrated contents (see above); the 
contribution of the micro-phytoplankton is extremely low (< 20%) with rather featureless 
profiles. 

In eutrophic conditions (S8 and S9), a micro-phytoplankton maximum tends to 
develop near the surface, while the pico-C and nano-C profiles become uniform, with low 
relative values. The nano-phytopl ank ton is predo minan t everywhere in the water column for 
intermediate trophic regimes, (i.e., when C sur f is around 0.2-0.5 mg mf 3 ; see S5 and S6). 

Parameterizing the micro-C, nano-C, and pico-C profiles is made as above for C, 
according to equation (7). The profiles reconstructed with the appropriate parameters (Table 
5) are displayed in Figures 8d- f. They reproduce fairly well the profiles derived from the 
statistical analysis. 

The same profiles, restored in their actual concentration and depth values, are 
displayed in Figure 9. The absolute magnitude of the pico-phytoplanktonic peaks remains 
about the same (0.09-0.1 1 mg nf 3 ) for the first four classes (SI to S4), in contrast to the nano- 
phytoplanktonic peaks which steadily increase (0.08 to 0.14 mg m' 3 ). 

4.2.2 Micro-, nano-, and pico-phytoplankton in well-mixed waters 

Global ocean 

About one third of the mixed waters stations (184 over 479) were located in the 
Southern Ocean (south of 60°S). These data are examined separately for reasons explicated 
thereafter. The 295 remaining stations from the global ocean (Southern Ocean excluded) are 
analyzed according to the same trophic categories as defined above (Ml to M5), with respect 
to their vertically integrated content, and then with respect to the shape of their profiles. 

The integrated contents of micro-C, nano-C, and pico-C, expressed as percentage of the 
total content, are displayed in Figure 10a (to be compared with Figure 7a). The figure shows 
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there is no notable difference between the mixed and stratified waters regarding their 
respective proportions of pico-, nano-, and micro-phytoplankton contents, except a lower 
contribution of pico-phytoplankton in oligotrophic conditions (30% vs 50% in mixed and 
stratified waters, respectively). 

As expected from the profiles of total chlorophyll a (Figure 5b), the individual profiles 
for pico-, nano-, and micro-phytoplankton (Figures lla-c) are also rather featureless, 
essentially uniform within the euphotic layer, and decreasing below (C, from 1 to 2). 
Regarding the size composition, they bear some resemblance with those in stratified waters. 
Indeed, in low C waters (C sur f < 0.4 mg m' 3 ), the micro-phytoplankton fraction is at its 

minimum (10% of Czcu), while the pico- and nano-phytoplankton dominate (with a smooth 
peak near £ = 0.8). The latter is relatively more abundant (50% vs: 30%) in contrast to what 
happens in stratified waters. In eutrophic regimes, the micro-phytoplankton fraction largely 
dominates (more than 80%), and this dominance is similar to what was described for 
stratified systems (Figure 7a). Note that, in such high C stratified systems, the depth of the 
euphotic zone is reduced to such a pomt that waters could be equally considered as mixed 
ones (according to Z eu /Z m ). 

Particular case of the Southern Ocean 

Southern (mixed) waters are treated separately in this study. A preliminary analysis of 
the data set showed that, contrarily to temperate latitudes where a C increase was generally 
associated with a Fuco (diatoms) increase, in the Southern Ocean, 1 9'-HF (prymnesiophytes) 
increase was usually observed during bloom periods. About 80% of these data originate from 
the Ross Sea. This area is known to sustain Phaeocystis blooms in relatively mixed 
conditions, whereas diatom blooms occur in stratified areas \Goffart et ah, 2000; Boyd, 2002; 
Arrigo et al., 2003]. Studies conducted in other southern areas also reported that 
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phytoplankton blooms can be dominated either by nano-phytoplankton species (e.g., 
Phaeocystis spp., cryptophytes, and nanoflagellates) or by diatoms [ Moline and Prezelin, 
1996; Claustre et al., 1997]. More generally, these observations and the subsequent partition 
of the data set (mixed waters) into two subsets (global and southern waters) are consistent 
with a recent body of literature which emphasizes the particular bio-optical status of polar 
waters [ Stramska et al., 2003; Cota et al., 2003], possibly related to taxonomic composition 
differences [Sathyendranath et al., 2001], when compared to temperate or tropical 
environments. 

The partition of the data set between southern and temperate (i.e., global without 
southern) waters appears a posteriori supported by the examination of Figure 10b (to be 
compared with Figure 10a). Indeed, in southern waters, the relative abundance of the three 
size classes is significantly different, and the pattern of the nano- and micro-phytoplankton 
contributions is reversed. Actually, micro-phytoplankton dominate at low C SUI f and decrease 
to the profit of nano-phytoplankton when C SU rf increases. At C surf greater than 2 mg m' 3 , nano- 
phytoplankton represent up to 60 to 80% of the biomass in southern waters (while diatoms 
would represent about 90% of the biomass in temperate waters). The pico-phytoplankton 
biomass, which never exceeds 10% of the total chlorophyll a content, is distinctly lesser in 
southern waters than in any other areas. 

The vertical profiles in the Southern Ocean (Figures 1 1 d-f) present a relatively uniform 
distribution (down to 2 Z eu ) for the nano-phytoplankton assemblages, and also show that this 
population is preponderant. A micro-phytoplankton deep maximum seems to develop above 
the euphotic depth, in low Csurf waters (Ml class). A closer examination of the corresponding 
profile shows a relatively high variability around the mean maximum value 0.8 Cze„ ± 
0.4). 
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Parameterization in mixed waters 


In the mixed regime, the average proportion of each phytoplankton class within the 
euphotic layer is dependent on the geographical location (global vs southern) and on the 
trophic status, as a function of C sur f. Furthermore, with the exception of a few situations, the 
vertical distribution of the micro-C, nano-C, and pico-C generally present weak changes with 
depth, at least within the euphotic layer. It is reasonable, therefore, to use the results 
presented in Figure 10 as a basis for inferring, from the sole knowledge of C su rf, the vertical 
distribution of micro-C, nano-C, and pico-C with depth. The fractional C is simply calculated 
by multiplying C sur f by the corresponding proportion, and this concentration is thereafter 
assumed to be homogenous from the surface down to Z eu - The proportions for each mixed 
trophic category are tabulated in Table 6. 

4.2.3 Validation of the method for predicting the micro-C, nano-C and pico-C profiles 

Vertical distribution 

The fraction of the HPLC database (20%, including 483 micro-C, nano-C and pico-C 
vertical profiles), which was kept for the validation process, is presently used (Figure 2). The 
euphotic depth is computed from the C SU rf value, and compared to Z m , to determine if the 
water column is stratified or mixed. For stratified waters, the C sur f value is injected into the 
relevant parameterization to generate the dimensionless vertical profiles which are then 
restored in physical units (by using the actual Z eu and Cz= u values). For mixed waters, 
vertical profiles are simply derived by extrapolating to depth the product of Csurf by the 
relative proportion of each group. By using these profiles, the column-integrated content (for 
micro-, nano-, and pico-phytoplankton) can be predicted (for the 0-Z eu interval), and directly 
compared to the measured contents (Figure 12). For the micro-C and the nano-C (Figures 
12a-b), the points are rather regularly distributed along the 1:1 line over a large range of 
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concentrations 0.5-100 mg m~ 2 for micro-phytoplankton, and ~ 2-100 mg m 2 for nano- 
phytoplankton). The distribution is apparently less satisfying for pico-phytoplankton, 
probably because the biomass of this group varies over a much more restricted range (2-8 mg 
m' 2 predicted vs 1.5-15 mg m' 2 measured). 

For stratified waters, the parameterization is further evaluated through the analysis of 
two descriptors of the profile shape: the maximum value of the micro-C, nano-C and pico-C 
profiles, and the depth Z max of this maximum (Figures 12d-j). The same comments given for 
Figures 12a-c apply here. The agreement between measured and predicted shape parameters 
is satisfactory as far as the micro- and the nano-phytoplankton are concerned. The situation is 
less favorable for pico-phytoplankton, but this is largely because of a weaker 
parameterization of this class in southern waters. Were these waters disregarded, the quality 
of the parameterization would be nearly similar for the three groups. 

Surface layer 

This validation uses the separate databank (database #2, see Figure 2) which includes 
4238 surface samples. The actual C concentrations (sampled between 0 and 10 m) are used as 
C sur f. The respective proportions of micro-, nano-, and pico-phytoplankton are retrieved for 
each C sur f, by using the results of Figure 7c for stratified waters, and of Figure 10 for mixed 
waters. The micro-C, nano-C, and pico-C are then computed by multiplying the 
corresponding proportions with C slff f, and finally compared to the measured values (Figure 
13). For the global ocean (southern waters excluded), the agreement is excellent over a wide 
range of C sur f values for the micro- and nano-phytoplanktonic groups, and also for the pico- 
phytoplankton, despite a much restricted range of variation. Again, the parameterization is 
less satisfying for southern waters. 
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4.2.4. Example of application of the proposed parameterization 

The relevance of the parameterization developed as part of this study can be shown by 
mapping the distribution of the three phytoplankton size classes at the global scale. A 
SeaWiFS monthly composite (for June 2000) is taken as an example. The C value associated 
with each pixel (Csurf) is used to infer the proportions, as well as the vertical distribution, of 
the biomass corresponding to the three phytoplankton size classes. A preliminary partition of 
the ocean into mixed and stratified waters is necessary before using the appropriate 
parameterization (Tables 5 and 6). The results display the typical phytoplankton signatures of 
the main oceanic entities (Figures 14a-c). 

Pico-phytoplankton appear to be the dominant group within the subtropical gyres 
where they can represent up to 50% of the biomass. In high latitude environments, however, 
their contribution is less than 10%, especially in the northern hemisphere. The reverse 
situation is observed for micro-phytoplankton, which always represent less than 20% 
(sometimes less than 10%) in the gyres, and, in contrast, 40-70% of the biomass in the North 
Atlantic (latitudes greater than 40°N), where they form the dominant size class at the end of 
the spring bloom. In between, nano-phytoplankton represent an extremely stable biomass 
everywhere (30%-50% of the total). 

The analysis can be further extended by considering the vertically integrated content 
of micro-C, nano-C and pico-C within the euphotic layer (Figures 14d-f). Micro- 
phytoplankton (likely diatoms) definitely present the largest range of biomass variation with 
values close to 0 in the center of the subtropical gyres and up to 80 mg m" 2 in the North 
Atlantic Ocean. Pico-phytoplankton biomass varies within a narrow range of concentrations 
(7-8 mg m' 2 ). As a result of its relative stability, the nano-phytoplankton distribution typically 
mimics that of the total biomass. 
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5. Conclusions 


Although the MB89 study was based on a substantial amount of data (3806 profiles), 
thus leading to rather robust results, it is still appropriate to test this robustness by considering 
an independent data set. Such a test is particularly required if an extension of the same 
approach is envisioned, dealing with other pigments, or groups of pigments, instead of the 
total chlorophyll a alone. 

The results of this test are highly satisfactory. The relationships between C sur f and the 
total chlorophyll a column-integrated content in both stratified and mixed waters of this study 
coincide perfectly with those obtained in MB89. This comforting agreement calls for two 
comments. The methodological change when determining C has no discernible impact upon 
the results, which is an indirect consolidation of the previous measurements, which were 
mainly based on the fluorometric method. Such global relationships testify to the existence of 
robust global trends over a wide range of trophic conditions. The scattering of the points (see 
Figure 3), however, simultaneously demonstrates that any relationship of that sort is not as 
well defined when a restricted range of C is considered [e.g. discussion in Siegel et ai, 2001]. 

Regarding the regular and progressive evolution of the shapes of the vertical profile 
according to C sur f, the strong resemblance (actually a close numerical agreement) between the 
previous results and those presently obtained is very encouraging. This confirmation 
reinforces the notion of a continuity in the relationship between C sur f and the featured vertical 
distribution typical of stratified waters (replaced by uniform distributions in mixed waters). 
Such a continuity can be clearly discerned when the statistical analysis is performed by using 
dimensionless quantities, for both the depth and the concentration, whereas it remains 
obscure when the physical quantities are directly considered. For modeling purposes, 
including those for primary production based on satellite surface chlorophyll data, such a 
continuity is helpful and easy to handle horn a numerical point of view. 
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The realization of the continuity is also supported by the evidence of regular changes 
within the pigment composition, the continuum of algal cell size, and the vertical layering of 
the dominant populations. The results dealing with the phytoplanktonic assemblages and the 
predominance of algal groups according to the trophic conditions, confirms existing 
qualitative descriptions, generally based on local observations. Though, owing to the 
statistical analysis, the present study extends the previous knowledge, and permits its 
generalization to most oceanic situations. 

The proposed statistical relationships, as all statistical products, may fail on a case-by- 
case basis, and their ubiquity is certainly doubtful in high latitudes. Indeed, polar zones are 
largely underrepresented in the present data set, with the exception of the Ross Sea which has 
received a particular attention because of its relatively high productivity. The hydrological 
and climatological conditions, prevailing in southern waters, have been found to favor 
particular phytoplankton assemblages with respect to those growing in temperate regions. A 
special statistical analysis performed on a more specific data set must, therefore, be 
considered for these particular waters, as well as for Arctic waters not represented here. 

Finally, this study constitutes a first step in the synthesis of HPLC data collected 
during the last decade. These data, originating from many sources, provide an initial 
framework to comprehensively describe the distribution of the phytoplankton communities at 
global scales, which has become increasingly needed to appraise how the algal biomass 
interacts with carbon fluxes. The generalized trends in the distribution of phytoplankton size 
classes according to the trophic status are not only confirmed in this study at global scales, 
but are also parameterized in a practical way; it thus becomes possible to propose a first 
assessment of the phytoplankton community composition for the world ocean from ocean 
color satellite observation (Figure 14). The transformation of “maps of phytoplankton size 
class biomass” into “maps of primary production specific to each size class”, through the use 
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of bio-optical models, appears now feasible (to the extent that the photosynthetic parameters 
typical of the different phytoplankton size classes can be selected). Such maps of “specific 
primary production” could in turn be very useful to improve the parameterization and the 
validation of recent biogeochemical models which include different phytoplankton 
compartments [ e.g . Moore et al., 2002; Aumont et al, 2003]). 
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Figures captions 


Figure 1: Geographic location of the 2419 stations used in the present study. The size of the 
black dots indicates the number of stations within a square of 10° longitude by 10° latitude. 
The small dot corresponds to less than 10 stations; the medium to 10-100 stations; and the 
large to more than 100 stations. 

Figure 2: Flowchart detailing the partition and the use of the two data sets for data analysis as 
well as for validation of the developed parameterization. 

Figure 3: Relationship between the average total chlorophyll a concentration within the 
surface layer, Csurf, and the total chlorophyll a integrated content within the euphotic 
layer, {C )^ , for stratified (a) and mixed (b) waters. The solid lines represent equations (8) 
while thfe dashed lines correspond to the regression lines models of MB89. 

Figure 4: C vertical profiles for stratified waters. Average dimensionless C profiles for each 
trophic category (a), modeled dimensionless C profiles for each trophic category (b), rescaled 
C profiles for categories SI to S4 (c), and rescaled C profiles for categories S5 to S9 (d). For 
rescaled profiles of each trophic category, the position of Z eu is identified by an arrow 
associated with the symbol relevant to the trophic category. 

Figure 5: C vertical profiles for mixed waters. Average dimensionless C profiles for each 
trophic category (a), rescaled C profiles for categories Ml to M3 (b), and rescaled C profiles 
for categories M4 and M5 (c). For rescaled profiles of each trophic category, the position of 
Z cu is identified by an arrow associated with the symbol relevant to the trophic category. 

Figure 6: Comparison of parameterizations developed in the present study and in MB89 to 
retrieve the C vertical profile from C sur f. Three criteria are compared. Integrated content 
within the euphotic layer (c) (a), maximal concentration Cmax (b) and depth of the 

maximal concentration Z max (c). 

Figure 7: Change in the size structure of the phytoplankton community as a function of C SU rf 
for stratified waters. Average proportion (%) of micro-C, nano-C and pico-C, within the 
euphotic layer (a). Average integrated content of micro-C, nano-C, pico-C and C within the 
euphotic layer (b). Average proportion (%) of micro-C, nano-C and pico-C, within the surface 
layer (c). The vertical bars represent ± 1 sd around mean values. 

Figure 8: Average (a-c) and modeled (d-f) dimensionless vertical profiles of rmcro-C, nano-C 
and pico-C in stratified waters. 

Figure 9: Average rescaled vertical profiles of micro-C (a-b), nano-C (c-d) and pico-C (e-f) 
in stratified waters. For each trophic category, the position of Z cu is identified by an arrow 
associated with the symbol relevant to the trophic category. 

Figure 10: Change in the size structure of the phytoplankton community as a function of C SUI f 
for mixed waters. Average proportion (%) of micro-C, nano-C and pico-C, within the 
euphotic layer for global mixed waters, southern waters excluded (a) and for southern mixed 
waters (b). The vertical bars represent ± 1 sd around mean values. 



Figure 11: Average dimensionless vertical profiles of micro-C, nano-C, and pico-C in mixed 
waters. Data are split into global mixed waters (excluding southern waters) (a-c), and southern 
mixed waters (d-f). 

Figure 12: Comparison of measured and computed characteristics of vertical profiles of 
micro-C, nano-C and pico-C. Integrated content within the euphotic zone for stratified and 
mixed waters (a-c); maximal concentration for stratified waters (d-f); depth of the maximal 
concentration for stratified waters (g-i). 

Figure 13: Comparison of measured and computed micro-C (a), nano-C (b) and pico-C (c) 
for surface waters. 

Figure 14: Global distribution of the size structure of the phytoplankton community for June 
2000 (SeaWiFS composite). Proportion (%) of micro-C, nano-C and pico-C (a-c), and 
integrated content within the euphotic layer (mg m' 2 ) of micro-C, nano-C and pico-C (d-f). 


Table 1: Symbols used in the present study and their significance. 


Symbol significance 


£ Total chlorophyll a (divinyi-chlorophyll a + monovinyl-chlorophyll a + 

chlorophyllid a ) concentration, mg m' 3 
PAR Photosynthetically available radiation, W m' 2 

K Attenuation coefficient for downwelling PAR irradiance, m' 1 

Zpd Depth of the surface layer, or penetration depth (= 1 / K), m 

Csurf Average integrated C within the layer O-Zpd, mg m 3 

z Geometrical depth, m 

„ Depth of the euphotic layer, defined as the depth where the PAR is reduced to 1% of 

its surface value, m 

£ Depth normalized with respect to Z eu : £ = z / Z cu , dimensionless 

(^)zcu Column-integrated content of pigment P within the euphotic layer, mg m 2 

Czeu Average column-integrated content of C within the euphotic layer, mg m 3 

P(Q Concentration of pigment P, at (the dimensionless depth) C mg m 3 

^ Concentration of pigment P normalized with respect to Czcu , p(Q = P(C) / Czcu , 

dimensionless 

micro-C Total chlorophyll a concentration associated to micro-phytoplankton, mg m 

nano-C Total chlorophyll a concentration associated to nano-phytoplankton, mg m" 3 

pico-C Total chlorophyll a concentration associated to pico-phytoplankton, mg m' 3 

[%micro] Proportion of biomass associated to micro-phytoplankton (%) 

[%nano] Proportion of biomass associated to nano-phytoplankton (%) 

[%pico] Proportion of biomass associated to pico-phytoplankton (%) 

Z m Depth of the mixed layer, m 

Pmax Maximum concentration on a considered P vertical profile, mg m" 3 

Zmax Depth Of P m ax 5 Ul 


Table 2: Biomarker pigments used in the present study and their taxonomic significance, associated size classe, and corresponding C to P ratio. 



Table 3: Statistical relationships between the average total chlorophyll a concentration 
within the surface layer, C SU rf (mg m' J ), and the total chlorophyll a integrated content within 
the euphotic layer, (C)^ (mg m' 2 ). The relationships obtained by MB89 are reproduced in 
italic. 

Stratified waters (n = 1821 profiles) 

For C torf < 1 mg m' 3 

(C} 2cu = 36.1 C sui( ° 357 (8a) (r 2 = 0.73, p< 0.001) (c)^ = 38.0 Cj 425 

For Cjurf > 1 mg m' 3 

( C )zcu = 37 ' 7 C surf°' 6li (8b) (r 2 = 0.70, p< 0.001) {C) Zai =40.2 Cj m 


Mixed waters (n = 598 profiles) 


( C )z.u =421C surf 0S38 (8c) (r 2 = 0.95, p< 0.001) (c)^^ 38.0 C 


sutf 


0.538 





Table 4: Trophic categories defined with respect to the average total chlorophyll a concentration within the surface layer, C sur f, and the associated relevant 
parameters. The parameters are derived from the calculations performed on the complete database #1, and are presented as averages and standard deviations 
(the latter are shown in parentheses). 

Stratified waters Mixed waters 




Table 5: Values of the 5 parameters of equation (7), obtained for the average dimensionless vertical profiles of C, micro-C, nano-C, and pico-C, 
for each trophic class of the stratified waters (S 1 to S9). 
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Table 6: Average proportions (%) of micro-C, nano-C, and pico-C within the euphotic layer, for the average C SU rf values of each trophic 
category for mixed waters (Ml to M5). 
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Figure 1: Geographic location of the 2419 stations used in the present study. The size 
of the black dots indicates the number of stations within a square of 10° longitude by 
10° latitude. The small dot corresponds to less than 10 stations; the medium to 10-100 
stations; and the large to more than 100 stations. 


Stratified waters 
(1821 profiles) 



Figure 2: Flowchart detailing the partition and the use of the two data sets for data analysis as well as for validation of the developed 
parameterization. 
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Figure 3: Relationship between the average total chlorophyll a concentration within 
the surface layer, C sur f, and the total chlorophyll a integrated content within the 
euphotic layer, (c) , for stratified (a) and mixed (b) waters. The solid lines represent 

equations (8) while the dashed lines correspond to the regression lines models of 
MB89. 
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Figure 4: C vertical profiles for stratified waters. Average dimensionless C profiles 
for each trophic category (a), modeled dimensionless C profiles for each trophic 
category (b), rescaled C profiles for categories SI to S4 (c), and rescaled C profiles 
for categories S5 to S9 (d). For rescaled profiles of each trophic category, the position 
of Z eu is identified by an arrow associated with the symbol relevant to the trophic 
category. 
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Figure 5: C vertical profiles for mixed waters. Average dimensionless C profiles for 
each trophic category (a), rescaled C profiles for categories Ml to M3 (b), and 
rescaled C profiles for categories M4 and M5 (c). For rescaled profiles of each trophic 
category, the position of Z eu is identified by an arrow associated with the symbol 
relevant to the trophic category. 
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Figure 6: Comparison of parameterizations developed in the present study and in 
MB89 to retrieve the C vertical profile from C sur f. Three criteria are compared. 
Integrated content within the euphotic layer (C) Z! . u (a), maximal concentration Cmax 
(b) and depth of the maximal concentration Zmax (c). 
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Figure 7: Change in the size structure of the phytoplankton community as a function 
of Csurf for stratified waters. Average proportion (%) of micro-C, nano-C and pico-C, 
within the euphotic layer (a). Average integrated content of micro-C, nano-C, pico-C 
and C within the euphotic layer (b). Average proportion (%) of micro-C, nano-C and 
pico-C, within the surface layer (c). The vertical bars represent ± 1 sd around mean 
values. 
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Figure 9: Average rescaled vertical profiles of micro-C (a-b), nano-C (c-d) and pico- 
C (e-f) in stratified waters. For each trophic category, the position of Z eu is identified 
by an arrow associated with the symbol relevant to the trophic category. 






Figure 10: Change in the size structure of the phytoplankton community as a function 
of Csurf for mixed waters. Average proportion (%) of micro-C, nano-C and pico-C, 
within the euphotic layer for global mixed waters, southern waters excluded (a) and 
for southern mixed waters (b). The vertical bars represent ± 1 sd around mean values. 
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Figure 11: Average dimensionless vertical profiles of micro-C, nano-C, and pico-C in 
mixed waters. Data are split mto global mixed waters (excluding southern waters) (a- 
c), and southern mixed waters (d-f). 
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Figure 12: Comparison of measured and computed characteristics of vertical profiles of micro-C, nano-C 
and pico-C. Integrated content within the euphotic zone for stratified and mixed waters (a-c); maximal 
concentration for stratified waters (d-f); depth of the maximal concentration for stratified waters (g-i). 
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Figure 13: Comparison of measured and computed micro-C (a), nano-C (b) and pico- 
C (c) for surface waters. 








Electronic supplement 

Appendix A: Summary of the data collection activities for both databases used in the present study. 

Al: Origin of the data included in database #1 




BENCAL 2002 Benguela Upwelling J. Ras 

BOFSD1183 1989 NE Atlantic http://www.bodc.ac.uk 

BOFSCD46 1990 NE Atlantic R. Barlow / http://www.bodc.ac.uk 

B0FSD1191 1990 NE Atlantic 

BOFS STERNA JCR2 1992 Bellingshausen Sea R. Barlow / http://www.bodc.ac.uk 
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METEOR 1995 Mediterranean Sea F. Vidussi 

MINOS 1996 Mediterranean Sea 

MIRC980CT 1998 Off Chile coast V. Stuart and E. Head 

MLML 91 1991 Off Iceland coast R. Bidigare 
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Appendix B: Average dimensionless profiles of C, micro-C, nano-C and pico-C, 

and their standard deviations 

Bl: Average dim ensionless C profiles for each trophic category (SI to S9) of 
stratified waters. The dashed lines represent the mean values ± sd. 
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B4: Average dim ensionless nano-C profiles for each trophic category (SI to S9) of 
stratified waters. The dashed lines represent the mean values ± sd. 


















Global mixed waters 


B6: Average dimensionless micro-C profiles for each trophic category (Ml to M5) of global (a-e) and 
southern (f-j) mixed waters. The dashed lines represent the mean values ± sd. 
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B7: Average dimensionless nano-C profiles for each trophic category (Ml to M5) of global (a-e) and 
southern (f-j) mixed waters. The dashed lines represent the mean values ± sd. 




B8: Average dimensionless pico-C profiles fqr each trophic category (Ml to M5) of global (a-e) and 
southern (f-j) mixed waters. The dashed lines represent the mean values ± sd. 









